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A pendulum prepared perfectly inverted and 
motionless is a prototype of unstable equilib- 
ria and corresponds to an unstable hyperbolic 
fixed point in the dynamical phase space. Un- 
stable fixed points are central to understand- 
ing Hamiltonian chaos in classical systems [1]. 
In many-body quantum systems, mean-field ap- 
proximations fail in the vicinity of unstable fixed 
points and lead to dynamics driven by quantum 
fluctuations [2, 3]. Here, we measure the non- 
equilibrium dynamics of a many-body quantum 
pendulum initialized to a hyperbolic fixed point 
of the phase space. The experiment uses a spin- 
1 Bose condensate [4—6], which exhibits Joseph- 
son dynamics in the spin populations that cor- 
respond in the mean-field limit to motion of a 
non-rigid mechanical pendulum [7, 8]. The con- 
densate is initialized to a minimum uncertainty 
spin state, and quantum fluctuations lead to non- 
linear spin evolution along a separatrix and non- 
Gaussian probability distributions that are mea- 
sured to be in good agreement with exact quan- 
tum calculations up to 0.25 s. At longer times, 
atomic loss due to the finite lifetime of the con- 
densate leads to larger spin oscillation amplitudes 
compared to no loss case as orbits depart from 
the separatrix. This demonstrates how decoher- 
ence of a many-body system can result in more 
apparent coherent behaviour. This experiment 
provides new avenues for studying macroscopic 
spin systems in the quantum limit and for inves- 
tigations of important topics in non-equilibrium 
quantum dynamics [9]. 

A pendulum initialized to a hyperbolic fixed point is 
metastable in the classical limit. Phase orbits passing 
close to these points have exponentially diverging peri- 
ods, and the orbits passing exactly through these points 
form a separatrix between librational and rotational mo- 
tion of the pendulum with an infinite period. If the pen- 
dulum is prepared perfectly in this orientation, the classi- 
cal equations of motion predict that it will not evolve. In 
reality, even if perfect preparation was possible, thermal 
fluctuations of the pendulum would perturb the pendu- 
lum from the metastable orientation and lead to oscil- 
lation. Even at zero temperature, unavoidable quantum 
fluctuations would lead to evolution [10, 11]. Although 
mechanical pendulums operating at the quantum limit 
are currently unavailable in the lab, it is possible to study 
quantum many-body systems that have similar dynami- 
cal behavior [12-14]. 

The focus of this work is spin-1 atomic Bose conden- 



sates [4-6] with ferromagnetic interactions tightly con- 
fined in optical traps such that spin domain formation 
is energetically suppressed. In this case, the non-trivial 
dynamical evolution of the system occurs only in the in- 
ternal spin variables, and the mean-field dynamics of the 
system can be described by a non-rigid pendulum similar 
to the two site Bose-Hubbard model [7, 8]. The system 
is fully integrable in both the quantum [15] and classi- 
cal [8, 16] limits, and exhibits a rich array of non-linear 
phenomena including Hamiltonian monodromy [17]. Fur- 
thermore, the condensate features a tunable Hamiltonian 
with a quantum phase transition that permits quench- 
ing of the condensate to highly-excited spin states. To- 
gether, these provide unique capabilities to explore non- 
equilibrium quantum dynamics that are not captured by 
mean-field approaches and can be solved exactly with 
Schrodinger's equation. 

In these experiments, we study the evolution of a 
quenched spin-1 condensate prepared in a metastable 
state corresponding to a hyperbolic fixed point in the 
spin-nematic phase space that ultimately evolves far be- 
yond the perturbative limit. The quantum solution of 
the problem at zero magnetic field yields intricate spin- 
mixing dynamics that exhibit non-linear quantum re- 
vivals [15] and a quantum carpet of highly non-Gaussian 
fluctuations [18]. At finite fields, the dynamics are simi- 
lar [19, 20], although they occur on a time-scale favorable 
for experimental observation. In both cases, the evolu- 
tion occurs along a separatrix of the phase space and is 
driven by quantum fluctuations that are absent from the 
mean-field theory solutions [16, 20]. 

The equilibrium states, domain formation and spin dy- 
namics of spinor condensates have been studied in many 
experiments [6, 13, 21-31]. In particular, observation of 
coherent spin oscillations have confirmed the mean-field 
pendulum model for small condensates [13, 23, 24]. Spin 
evolution has been previously observed from metastable 
spin states in many experiments [6, 22, 27-31], however, 
the experiments have not yet demonstrated spin dynam- 
ics in agreement with quantum calculations, except in 
the perturbative, low-depletion limit at very short times 
(where a Bogoliubov expansion around the mean field can 
be used) [26, 27, 31] or for conditions where the mean- 
field approach suffices. Here, by using low- noise atom 
detection techniques and careful state preparation, we 
are able to observe quantum spin dynamics that agree 
well with quantum calculations and demonstrate a rich 
array of non- Gaussian fluctuations. 

We begin by discussing the exact quantum model for 
spin-1 condensate small enough to be described by a sin- 
gle domain. The quantum states of the system can be de- 



2 




FIG. 1: Phase Space The spin state immediately following 
the quench is depicted on two relevant phase spaces of the 
spin-1 system: the {po,0 s } phase space (bottom-left) and the 
{S x ,Qyz,Qzz} spin-nematic Bloch sphere (right). A zoom- in 
of the hyperbolic fixed point at the pole is shown (top-left) 
with arrows indicating the orbit directions. The po,6 s phase 
space represents a Mercator projection of the {S x , Q yz , Qzz} 
sub-space. 

scribed in a Fock basis, |iVi, TVo, ^V-i), where Ni are the 
number of atoms in the three spin-1 Zeeman states. The 
spin dynamics, including the effects of a magnetic field, 
are governed by the interaction Hamiltonian [15, 19, 20]: 

H = A[(7V 1 -7V_ 1 ) 2 + (27V -1)(7V 1 +7V_ 1 ) 
+ 2a\a\ai(i-i + 2a|al 1 ao&o] 
+ q(N 1 +N_ 1 ). (1) 

Here, hi are the bosonic annihilation operators for the 
three spin states and Ni — a\ai. A and q(cx B 2 ) char- 
acterize the inter-spin and Zeeman energies, respectively. 
The spin-dependent binary collisions restrict the dynam- 
ical evolution to states that conserve both the total num- 
ber of atoms TV = Y^i Ni and the projection of angular 
momentum along the quantization axis M = N\ — N-±. 
Starting from the initial state |0, TV, 0), consisting of all TV 
atoms in the rrif = state, the evolution is constrained 
to final states of the form ^2 p c p \p, N — 2p,p). Hence, 
the solution to the quantum many-body problem is fully 
enumerated by the time-dependence of the Fock state 
amplitudes, c p (t). 

The semi-classical dynamics of the system take the 
form of a non-rigid pendulum [8]. Mean field states of a 
spin-1 condensates can be written as ip = (C+i, Co? C-i) T 
where Q = yj~pie %d % and pi = \d\ 2 = Ni/N are the frac- 
tional spin populations. The conservation of magneti- 
zation m = (TVi — N-i)/N constrains the populations 
p±i = (1 — po ± m)/2, and for the m = case that 
is relevant for these experiments, the spin dynamics are 



determined by the Hamiltonian: 

H = X'x 2 - A'(l - x 2 )cos6 s - qx (2) 

Here, x = (p - 1/2) /2 and S = +1 + 0_i - 20 are 
canonically conjugate variables and X' = 27V A. This 
Hamiltonian has the form of a classical non-rigid pen- 
dulum and is similar to the double-well Bose-Hubbard 
model that has been used to study Josephson effects in 
condensates. The Hamiltonian can also be written us- 
ing a phase space of the spin vector Si and nematic 
(quadrupole) tensor Qij matrix operators for the spin- 
1 system: H = X f ^ i S 2 + qQ zz /2. The phase spaces for 
both of these forms are shown in Fig. 1, where it is clear 
that the po,6 s phase space corresponds to a projection of 
the spin-nematic phase space. 

The initial state of the system following the quench, 
1 0,7V, 0), is indicated in the different phase spaces in 
Fig. 1 using quasi-probability distributions of the initial 
state determined from the quantum uncertainties [31]. In 
the spin-nematic space, the state corresponds to a mini- 
mum uncertainty state centered at the pole. The pole is 
a hyperbolic fixed point lying at the intersection of the 
separatrix that separates the librational and rotational 
orbits of the system. In the projected po> @s phase space, 
the distribution in po is tightly packed at the top of the 
phase space with random spinor phase. In the absence of 
quantum fluctuations, the state initialized the hyperbolic 
fixed point is non-evolving. However, quantum fluctua- 
tions populate a family of orbits that straddle the fixed 
point, and subsequent evolution leads to phase flow along 
the unstable manifolds of the separatrix. In the short 
term, this creates squeezed states with negligible change 
in po [31]. For longer times, the system evolves along the 
separatrix, which forms a closed homoclinic orbit in the 
spin-nematic space. 

We now turn to the experimental results. The exper- 
iment begins with a rubidium-87 condensate containing 
4 x 10 4 atoms, initialized in the / = 1, rrif = hyperfine 
state and held in a high magnetic field. The condensate is 
rapidly quenched by lowering the field, and the spin pop- 
ulations are measured for different evolution times. The 
experiment is repeated many times in order to acquire 
sufficient statistics to determine the full probability dis- 
tributions of the populations. The main results of the pa- 
per are shown in Fig. 2, which shows the measured proba- 
bility density of po = No/N versus evolution time, which 
is effectively a determination of the probabilities \c p \ 2 . 
The experimental results are compared with a quantum 
calculation using a spinor energy, 2A7V = — 2tt^ x 7.5 Hz, 
chosen to match the population dynamics. Both the ex- 
periment and quantum solutions exhibit population evo- 
lution that is in good overall agreement. In particular, 
both exhibit a long pause (80 ms) before any population 
evolution is apparent. After this pause, the spin popu- 
lation executes a regular damped oscillation. Population 
evolution from the metastable state is exponentially sen- 
sitive to initial population in the mp = ±1 states [27]. 
At the earliest evolution time studied (15 ms), the to- 



5 




time (ms) 

FIG. 2: Time evolution of spin populations. Probability 
density of the fractional population of the condensate in the 
ra/ = state, poj as a function of time. The curves show 
the mean, po (black line) and ± the standard deviation, a 
(blue lines), (a) Experimental data showing the results of 50 
runs at each evolution time placed into 40 bins. The mean 
and standard deviation curves have been smoothed using a 
cubic spline, (b) Quantum calculation using the initial atom 
number, magnetic field ramp, and atom loss rate measured 
in the experiment. The Fock state probabilities |c p | 2 , placed 
into 100 bins, are plotted. 



tal population in these states is measured to be < 30 
atoms which represents an upper bound limited by atom 
detection noise [31]. Initial populations at this level ef- 
fect the duration of the initial pause and first oscillation 
minimum, but not the overall character of the evolution 
[18] (see Supplemental Information). For evolution times 
beyond > 250 ms, it is necessary to include in the the- 
ory the effects of atomic loss due to the lifetime of the 
condensate r = 1.8 s, which is discussed in more detail 
below. 

It is clear that the mean and standard deviation are 
insufficient to fully characterize the distribution of po for 
both the experiment and theory, since for much of the 
evolution the mean does not pass through the highest 
probability density, and the asymmetry indicates a sig- 
nificant skew in the distribution. This point is reinforced 
in Fig. 3, which shows the full probability distributions 
for several evolution times, along with the theoretical 
predictions. The highly non-Gaussian nature of the dis- 
tributions provide compelling evidence of the quantum 
nature of the spin dynamics. The physical origin of these 
non- Gaussian fluctuations is dispersion of neighboring 
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FIG. 3: Full probability distributions of po. Evolution on 
the semi-classical phase space (left column) and histograms of 
the measured spin population, po (right) for different evolu- 
tion times after the quench, 15 ms, 130 ms, 140 ms, 170 ms, 
and 240 ms. a. The simulations use the semi-classical equa- 
tions of motion together with a quasi-probability distribution 
for the initial state. The mean value for po and 6 S are in- 
dicated with a black dot. b. The histogram bars for each 
evolution time depict the measured probability density of po 
for over 900 experimental runs, and the red line represents 
the simulation. The purple bar in each histogram represents 
the bin in which the mean of po is located. 
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orbits about the separatrix. Immediately following the 
quench, the distribution in po is tightly packed at the 
top of the phase space with random spinor phase. This 
state corresponds to a minimum uncertainty state of the 
spin-nematic subspace shown in Fig. 1 [31]. As evolution 
proceeds, the phase, # s , converges towards the separa- 
trix separating the librational and rotational trajecto- 
ries, and the population starts to evolve along it. The 
separatrix has a divergent period [8], and so the states 
disperse significantly due to the different evolution rates 
of nearby energy contours. It is this dispersion, together 
with the shape of the orbit, that gives rise to the highly 
non-Gaussian probability distributions. 

In order to characterize the non-Gaussian distribution, 
we determine several central moments, Uk = ((p — p) k ) 
from the data. The first six central moments are shown 
in Fig. 4 compared with the quantum simulation. Over- 
all the measured moments are in good agreement with 
the predicted moments from the simulation. The popu- 
lation revival in the second oscillation predicted from the 
simulation is clearly seen in the first four moments, but 
is less obvious in higher moments. 

We now turn to a discussion of the role of atomic loss 
in the dynamical evolution. The lifetime of the conden- 
sate r = 1.8 s is only a factor of 10 larger than the spin 
evolution timescale (^150 ms), hence one expects that 
loss plays an important role in the dynamics. We ex- 
plore this question in Fig. 5 where we compare quantum 
calculations without loss, quantum calculations includ- 
ing uncorrelated loss and the experimental data. Un- 
correlated atom loss is incorporated into the calculation 
using quantum Monte Carlo (QMC) techniques with the 
collapse operators d — y^T/rcti. The loss causes the 
overall magnetization M to execute a random walk with 
a restoring tendency towards M = and decreases the 
spinor dynamical rate, which scales as A oc TV -3 / 5 . (Sup- 
plementary Information) 

For the first 250 ms of evolution corresponding to the 
first spin oscillation, the effects of loss are not discernable 
between the two calculations, and the experimental data 
are in good agreement with both. Beyond 250 ms, there 
are significant differences between the two quantum cal- 
culations. The spin population of the calculation without 
loss nearly returns to the initial value and then experi- 



ences a long pause followed by complex multi-frequency 
oscillations. The calculation with loss however exhibits 
steady oscillations with one dominant frequency and a 
slowly decreasing amplitude centered on the ground state 
populations. In the semi-classical picture, the apparent 
damping of the calculation without loss derives from the 
dispersion about the separatrix in Fig. 3. The effect of 
loss is to eventually move the orbits away from the sepa- 
ratrix, which turns off this dispersion and leads to more 
regular oscillations. 

While the inclusion of loss into the model makes a sig- 
nificant improvement in the agreement of long term dy- 
namics (> 250 ms) with the experimental results, it is 
clear that this simple loss model is inadequate to fully 
replicate the measurements at longer time scales. While 
the experimental data and the simulations with loss are 
qualitatively similar, there is clearly more dissipation in 
the experiment as the amplitude of the oscillations damp 
more quickly and the standard deviation decreases. In fu- 
ture work, we intend to further investigate the damping 
of the spin dynamics and its connection to thermaliza- 
tion of isolated quantum systems subject to loss. Similar 
investigations are on-going using 1-D condensate systems 
[32-35] , and it will be interesting to explore the similari- 
ties and differences in these completely different systems. 
Finally, we believe that our results point the way to a host 
of fascinating explorations of out-of-equilibrium quantum 
spin systems [3, 9]. 

Methods 

We prepare a condensate of TV = 38, 500 ± 500 87 Rb 
atoms in the \f — l,ra/ = 0) hyperfine state in a high 
magnetic field (2 G). The condensate is tightly con- 
fined in an optical dipole trap with trap frequencies of 
250 Hz. To initiate dynamical evolution, the condensate 
is quenched below the quantum critical point by lowering 
the magnetic field to a value 210 mG and then allowed to 
freely evolve for a set time. The trap is then turned off 
and a Stern-Gerlach field is applied to separate the rrif 
components during 22 ms time- of- flight expansion. The 
atoms are probed for 400 /is with three pairs of orthog- 
onal laser beams, and the resulting fluorescence signal is 
collected by a CCD camera with > 90% quantum effi- 
ciency. 
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FIG. 4: Evolution of po central moments. The mean value of po, the standard deviation a, and the third through sixth 
central moments, Uk — {(p — p) k ), k — 3 — 6. The odd moments are on top, and the even moments are on the bottom. In 
each plot, the black markers represent the results of 50 experimental runs, and the error bars are estimated using a bootstrap 
method for the third moment and higher. The red curves are the prediction of the quantum simulation including loss. 
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FIG. 5: Long-term evolution of po. The top graph shows 
the quantum calculation without loss. The middle graphs 
shows the calculation including the effects of uncorrelated 
atom loss. The bottom graph shows the measured data. In 
each plot, the mean value is shown as a solid line, and the 
shaded envelopes correspond to the standard deviation 
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In this Supplementary Information, we provide an 
overview of the calculation methods used to compare the 
experimental results to theory and discuss the effect of 
spin impurities in the initial state. 



is estimated from population dynamics and then scaled 
according to the condensate lifetime. This simple model 
captures much of the early dynamics, but does not pro- 
duce as much damping as observed experimentally. 



I. SECOND QUANTIZED CALCULATIONS 

A. Exact quantum calculations 

The first calculation method is the exact second quan- 
tized form of the Hamiltonian. The second quan- 
tized form operates on a Fock basis represented by 
|7Vi, 7V , N-i) where Ni is the number of particles in 
the rrif state given in the index. The spin mixing 
Hamiltonian conserves both the total atom number N = 
TVi + N + 7V_i and the magnetization M = N x - N_ x . 
The Hamiltonian in this basis is given by: 

Hsma = A[(iVx - iV_!) 2 + (2N - l)(iVi + N-i) 
+ 2a\a\aia-i + 2ajal 1 ao&o] 
+ q(N 1 +N_ 1 ). (1) 

Similar to Ref. [1-5], we write the Hamiltonian in a Fock 
basis of A/", M, and /c, where No = N — 2k. k represents 
the number of pairs of atoms that have evolved from the 
rrif = state via spin-mixing. In matrix form operating 
on the vector of k coefficients this can be written: 

Hk,k' = 

{AM 2 + 2Xk(2(N - 2k) - 1) + q(2k + \M\)}5 kik > 
+ 2\{(k' + l)y/(N - 2k')(N - 2k' - 1)4,^+1 
+ k'y/{N - 2k' + 1)(N - 2k' + 2)4,fc'-i}. (2) 

where A = f / dr |0(r)| 4 , q = fi 2 B B 2 J (h 2 E HFS ), c 2 is 
the spin-dependent collisional interaction strength, and 
Ehfs is the ground state hyperfine splitting. A and q 
characterize the inter-spin and Zeeman energies respec- 
tively. This produces a tridiagonal matrix which can be 
exactly solved or numerically integrated. In our simu- 
lation we numerically integrate since the magnetic field 
and spinor energy hence q and A vary in time. The value 
for q(t) is modeled from experimental measurements of 
the magnetic field using microwave spectroscopy. A varies 
with the atom number as Af _3//5 .[18] The simplest way 
to account for this in the dynamical simulations is to 
make A a function of time. Since N decays approxi- 
mately exponentially with atom loss, the value for X(t) 



B. Quantum Monte Carlo 

A more rigorous calculation of the effects of atom loss 
is to use a quantum Monte Carlo simulation. The quan- 
tum Monte Carlo is implemented similar to Refs. [6-8]. 
The atoms are assumed to be lost one at a time and 
the process of losing an atom effectively measures its m f 
state and so the collapse operators are simply related to 
the annihilation operators for the modes of the conden- 
sate. The numerical integration of the k coefficients is 
performed with an effective Hamiltonian, 

H eff = Hsma - % - ^ C\C i 

i 

= Hsma - (3) 
It 

where Ci = ^/l/rdi are the collapse operators for each 
mode (i = —1, 0, 1) and r is the condensate lifetime. Dur- 
ing the time interval At of the numerical integration each 
atom has a probability e~ At ^ T of remaining. For each 
atom a random number in the range to 1 is generated 
to stochastically determine how many atoms to annihi- 
late in each mode. If this number is greater than e -A£ / r , 
then the appropriate collapse operator is applied to the 
state vector. The number of atoms for each mode is 
given by (Ni). After the collapse operators have been 
applied the k coefficients are renormalized and the next 
step of the numerical integration is performed with up- 
dated values for N and M. Results are obtained from the 
quantum Monte Carlo simulation by taking the average 
of quantum expectation values from many runs with the 
same initial conditions but a uniquely seeded sequence of 
random numbers to determine the annihilation probabil- 
ities. In effect, the results of the quantum Monte Carlo 
simulation are the average of many quantum trajectories. 

At first glance the quantum Monte Carlo is a daunting 
task since in general it should be necessary to use a ba- 
sis spanning every possible value of N(t), M(t), and k(t) 
which scales as A^ 3 . However the action of the collapse 
operators shifts the state vector from N(t) and M(t) to 
N(t + At) and M(t + At) while modifying the k coeffi- 
cients in a well characterized way. At any given step of 
the calculation there is only one value for A^ and M. So 
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for any step of the calculation the basis is proportional to 
N(t) and is completely described with the current values 
of TV, M, and the complex coefficients for the k index. 



II. MEAN FIELD CALCULATIONS 
A. Semi-classical simulations 

To make connections with the mean-field theory, we 
use a semi-classical technique together with quasi-proba- 
bility distributions (QPD) to regain the quantum statis- 
tics. The mean-field order parameter is represented by 
a complex vector (Ci,(o,C-i) T wnere & represents the 
amplitude and phase of the classical field for the mode 
associated with the m f state given in the index. A mean 
field analysis of the spin Hamiltonian [9-11] produces the 
dynamical equations [12]: 

= E 1 C 1 + c[( Pl + po - p_i)Ci + CoCi] (4a) 
= E C + c[( Pl + p_i)Co + 2CiC-iCo*] (4b) 
= £7.^-1+ c[(p_i+po-pi)C-i+CoCr] (4c) 

= 2AA, pi = Ni/N = |0 1 2 is the fractional pop- 
ulation of the rrif = i component, and E{ is the magnetic 
field energy for each mode. 

A convenient parameterization of the order parameter 
is given by 



where c 





Po + m 
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Ci 

Co 

C-i 

(5) 

where m = 0\ — 0-\ is the magnetization or Larmor 
precession phase, S = 0\ +9-i — 20o is the spinor phase, 
and m = (Ni — N_i)/N is the fractional magnetization. 

Using this parameterization, the dynamical equations 
reduce to [13]: 

2 C 

Po = -jT Po Vi 1 ~ Po) 2 - ™ 2 sin 6 S 

2q 2c 

9 ' = ~n + n x 

(1 - 2p ) + = — cos a 

V(l - po) 2 - m 2 



(6) 



po and S are canonically conjugate variables and the 
dynamical phase space is defined by the corresponding 
spin energy functional of the condensate: 

£ ^m 2 +cp (1 - po) + - po) 2 - m 2 cos S 

+pra + g(l -p ). (7) 

m and m are also canonically conjugate variables. How- 
ever, m is cyclic because it does not appear in the energy 
functional, and thus m is conserved. 



1. SU(3) 

A spin-1 condensate can be described using the SU(3) 
symmetry group with the Cartesian dipole-quadrupole 
decomposition of the Lie algebra su(3) as a basis. There 
are three dipole (or angular momentum) operators S a , 
and nine quadrupole operators Q a b which are moments 
of the quadrupole tensor ({a, b} G {x,y,z}). Only five 
of quadrupole operators are linearly independent since 
the quadrupole tensor is symmetric and traceless. These 
operators are defined as [14, 15]: 



S a = -ihcabcbtb, 
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ab 



-b%-b\b a 



b u c 
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(8) 



(9) 



summing over repeated indices where 

bt = (-a^+al^A/2 
b\ = i(al ± +0,1^/^2 
bt = at 



In Figure 3 of the main text we depict an SU(2) sub- 
space of this SU(3) as {S x , Q yz ^Qzz}- Just as in [16], the 
exact subspace is {S x ,Q yz , (Q zz — Q yy )}- However, our 
initial state of m/ = makes the off-diagonal elements 
(ex ^/NiN-i) of (Q zz — Q yy ) small and always aver- 
age to zero over a Larmor precession cycle, so ignoring 
these terms (Q zz — Q yy ) ~ \Qzz- Additionally, Q zz is 
linear with A^o ( Q Z z = |(^ — 3iVo)), making for an easy 
mapping to po = Nq/N. 

For the initial state of all atoms in the rrif = level, 
1 0,7V, 0), we have the following expectation values for 
the variances: ((A^) 2 ) = ((A^) 2 ) = ((AQ^) 2 ) = 
((AQ yz ) 2 )=N 



2. Quasiprobability distributions 

It is well known that the perfect rrif = state does 
not evolve from these mean-field equations. Some seed- 
ing is necessary. This seeding comes from the quantum 
noise of the initial Fock state, |0, A, 0). In order to pro- 
duce the necessary QPD, we note that the |0, A, 0) Fock 
state is a minimum uncertainty coherent state with S x , 
Q yzi S y , and Q xz all having an expectation value of zero 
and variance of N. Here, S a ({a} G {x,y,z}) are the 
spin-1 dipole (or angular momentum) operators, and Q a 6 
({a, b} G {x, y, z}) are moments of the spin-1 quadrupole 
(nematic) tensor. Generating a sample of random num- 
bers for each of these distributions, we can then find val- 
ues for the initial fi, Co, and £-1 (or po and S ) to plug 
into the spinor dynamical equations for the semi-classical 
simulation [16]. The conversion equations are: 

H~ Qyz 



tanx+ 



$>x H~ Qx 
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Spin evolution from the metastable state is a para- 
metric amplification process whose early-time dynamics 
are exponentially sensitive to initial population in the 
rriF = ±1 states [17]. Hence, any impurities in the initial 
state preparation will certainly effect the timescale of the 
initial pause and first oscillation minimum. Importantly 
though, as shown in [2] , the overall character of the evo- 
lution, including the intricate evolution of the quantum 
spin fluctuations, is robust to impurities even up to the 
few percent level, which is an order of magnitude larger 
than our measured bound. 
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FIG. 1: Analysis of impurities. Upper bound from direct 
measure of impurities at 15 ms of evolution, gray shaded re- 
gion. Measurement of first minimum of po, orange shaded 
area. Simulation with polluting atoms for a spinor dynamical 
rate of c/(2nh) = 5.0, 5.5, 6.0, 6.5, 7.0, 7.5, and 8.0 Hz (blue 
lines top to bottom). The values between 7-8 Hz are most 
consistent with the long term dynamics and are shaded blue. 
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m = pi- p-i 



where \± = #±i ~~ 




^X H~ QxZ 

cos x+ 



$x Qxz 
COSX- , 



Each of these sample wavefunctions is numerically in- 
tegrated using the dynamical equations to yield po (t) and 
3 (i), which are then plotted on the semi-classical phase 
space as in Fig. 2 in the main paper. 



III. EFFECT OF INITIAL STATE IMPURITIES 

In this section, we discuss the effects of spin impurities 
in the initial state preparation, and the impurity lim- 
its determined by the experiment. The experiment uses 
condensates containing 4 x 10 4 atoms, initialized in the 
/ = l,ra/ = hyperfine state with measured impurities 
in the m/ = ±1 states (limited by the atom counting 
noise [16]) below < 30 atoms or 0.1% of the total popu- 
lation. 



In order to analyze the quantitative effect of impuri- 
ties, we perform simulations with two types of impurities: 
an initial non-zero magnetization and an initial non-zero 
number of pairs of m/ = ±1 atoms. The results of these 
calculations are nearly identical for the same number of 
impurity atoms with the non-zero magnetization results 
shown in Figure 1 for various levels of impurities and for 
a range of spinor dynamical rates determined from the 
long time evolution of the experiment (blue shaded re- 
gion) as well as several other values (blue lines). These 
are compared to experimental measurements in order to 
ascertain an upper bound on the impurities in the exper- 
iment. The first time at which the atom populations are 
measured is 15 ms after the beginning of the magnetic 
quench, which provides an upper bound on the impuri- 
ties present at t = 0. The population in the = ±1 
states at this time is < 30 atoms, which is shown as the 
gray shaded region in Fig. 1. Also plotted in Fig. 1 is 
the measured time that the rrif = population reaches a 
minimum value. This time is exponentially sensitive to 
impurity atoms and, for both magnetization and pair im- 
purities, reduces similarly for the same number of impu- 
rity atoms. The experimental measurement of this time 
plus and minus one standard deviation is shown as the 
orange shaded region. The shaded regions overlap only in 
the limit of very little pollution of the initial state. The 
overlap region is consistent with no pollution and is in- 
consistent with pollution of the magnetization of greater 
than 5-10 atoms and pair pollution of greater than 3-5 
pairs of ±1 atoms, even for significantly different spinor 
dynamical rates than the evolution suggests. While it is 
conceivable to trade off between spinor dynamical rate 
and pollution, the dynamics of the quadrature squeez- 
ing measurement reported previously [16] indicates that 
the value of spinor dynamical rate estimated from the 
long term population dynamics is more consistent with 
the available data. The analysis presented here along 
with the non- Gaussian nature of the fluctuations and the 
squeezing dynamics reported previously make an effec- 
tive argument for the initial state preparation producing 
a very pure m/ = state. 
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